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Abstract. Our aim in this paper is to study the Cahn-HiUiard equation with singular 
potentials and dynamic boundary conditions. In particular, we prove, owing to proper 
approximations of the singular potential and a suitable notion of variational solutions, the 
existence and uniqueness of solutions. We also discuss the separation of the solutions 
from the singularities of the potential. Finally, we prove the existence of global and 
exponential attractors. 



plays an essential role in materials science as it describes important qualitative features 
of two-phase systems related with phase separation processes. This can be observed, e.g., 
when a binary alloy is cooled down sufficiently. One then observes a partial nucleation 
(i.e., the apparition of nucleides in the material) or a total nucleation, the so-called spin- 
odal decomposition: the material quickly becomes inhomogeneous, forming a fine-grained 
structure in which each of the two components appears more or less alternatively. In a 
second stage, which is called coarsening, occurs at a slower time scale and is less under- 
stood, these microstructures coarsen. We refer the reader to, e.g., [5], [6], [28], [29], [31] . 
[32] . [lO] and [H] for more details. Here, u is the order parameter (it corresponds to a 
(rescaled) density of atoms) and fj, is the chemical potential. Furthermore, / is a double- 
well potential whose wells correspond to the phases of the material. A thermodynamically 
relevant potential is the following logarithmic (singular) potential: 



although such a potential is very often approximated by regular ones (typically, /(s) = 
— s). Finally, n is the mobility and a is related to the surface tension at the interface. 
This system, endowed with Neumann boundary conditions for both u and (meaning 

that the interface is orthogonal to the boundary and that there is no mass fiux at the 
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1. Introduction 



The Cahn-HiUiard system 




dfU = nAr^fi, K > 0, 

fi = —aA^u + f{u), a > 0, 
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boundary) or with periodic boundary conditions, has been extensively studied and one 
now has a rather complete picture as far as the existence, uniqueness and regularity of 
solutions and the asymptotic behavior of the solutions are concerned. We refer the reader, 
among a vast literature, to, e.g., P, [10], [E], p], ^9\, [20], [26], [30], [3l], [37], [3H], 
[39], [10], [H], [16], [50] and [51]. 

Now, the question of how the process of phase separation (that is, the spinodal de- 
composition) is influenced by the presence of walls has gained much attention recently 
(see [21], [22], [27] and the references therein). This problem has mainly been studied for 
polymer mixtures (although it should also be important in other systems, such as binary 
metallic alloys): from a technological point of view, binary polymer mixtures are particu- 
larly interesting, since the occurring structures during the phase separation process may 
be frozen by a rapid quench into the glassy state; micro-structures at surfaces on very 
small length scales can be produced in this way. 

In that case, we again write that there is no mass flux at the boundary. Then, in order 
to obtain the second boundary condition, following the phenomenological derivation of 
the Cahn-Hilliard system, we consider, in addition to the usual Ginzburg-Landau free 
energy 

(1.3) ^GL(n,Vu)= [ i^\VM'' + F{u))dx, 

where F' = f and Q is the domain occupied by the material (the chemical potential fi 
is defined as a variational derivative of \E'gl with respect to u), and assuming that the 
interactions with the walls are short-ranged, a surface free energy of the form 

(1.4) <iIriu,Vru) = [ {^\Vru\^ + G{u))dS, ar > 

Jr 2 

(thus, \E' = "^GL + is the total free energy of the system), where F is the boundary 
of Q and Vr is the surface gradient. Writing finally that the system tends to minimize 
the excess surface energy, we end up with the following boundary condition: 

(1.5) -dtu — arA-pu + g(u) + adnU = 0, on F, 

a 

where Ar is the Laplace-Beltrami operator, (9„ is the normal derivative, g = G' and 
(i > is some relaxation parameter, which is usually referred to as dynamic boundary 
condition, in the sense that the kinetics, i.e., dtu, appears explicitly. Furthermore, in the 
original derivation, one has G{u) = ^a^u^ — bru, where ar > accounts for a modification 
of the effective interaction between the components at the walls and 6r characterizes the 
possible preferential attraction (or repulsion) of one of the components by the walls (when 
6r vanishes, there is no preferential attraction). We also refer the reader to [3] and [23j 
for other physical derivations of such dynamic boundary conditions, obtained by taking 
the continuum limit of lattice models within a direct mean-field approximation and by 
applying a density functional theory, to [H] for the derivation of dynamic boundary 
conditions in the context of two-phase fluids flows and to ^47j and [H] for an approach 
based on concentrated capacity. 
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The Cahn-Hilliard system, endowed with dynamic boundary conditions, has been stud- 
ied in P, [21], [22], [12], [IS] and [SUj for regular potentials / and g. In particular, one 
now has satisfactory results on the existence, uniqueness and regularity of solutions and 
on the asymptotic behavior of the solutions. 

The case of nonregular potentials and dynamic boundary conditions is essentially more 
complicated and less understood. Indeed, to the best of our knowledge, even the exis- 
tence of weak energy solutions has only recently been established in that case, under the 
additional restriction that the boundary nonlinearity g has the right sign at the singular 
points ±1, namely. 



(see [25]; see also [7J where sign conditions are considered in the context of the Caginalp 
phase-field system). Furthermore, the questions related with the longtime behavior of the 
solutions (e.g., in terms of global attractors or/and exponential attractors) have not been 
considered in the literature. 

The aim of the present paper is to give a thorough study of the singular Cahn-Hilliard 
problem endowed with dynamic boundary conditions. As we will see below, the main 
difficulty here lies in the fact that the combination of dynamics boundary conditions and 
of singular potentials can produce additional strong singularities on the corresponding 
solutions close to the boundary (especially in the case where the sign condition fll.6p is 
violated). In that case, even the simplest ID stationary problems may not have solutions 
in a usual (or distribution) sense (due to the jumps of the normal derivatives close to the 
boundary produced by the singularities, see Example 6.2). 

Nevertheless, we can construct a sequence of solutions of regular approximations of our 
singular problem which converges to a unique trajectory which is then naturally identified 
with the "solution" of the limit singular problem. As already pointed out, this trajectory 
may not be a solution of our equations in the usual (distribution) sense, so that the notion 
of a solution must be properly modified. To do so, we consider, in the spirit of [1] (see 
also [12j), the variational inequality associated with the problem and define a (variational) 
solution in terms of this variational inequality, see Section 3 for details. 

Of course, important questions are when the solution thus defined is a usual distribution 
solution of the equations and which additional regularity one can expect from such a 
variational solution. Actually, we prove that the variational solutions are always Holder 
continuous in space and are solutions in the usual sense if they do not reach the pure 
states on the boundary, namely, if 



for almost all {t, x) G x F. One possible condition which guarantees that condition 
(ll.7p holds is exactly the aforementioned sign condition (11.61) (see Proposition 4.5). Al- 
ternatively, this condition is always satisfied if the singularities of the nonlinearity / are 
strong enough, namely, if 



see Section 4. Furthermore, using some proper modification of the Moser iteration scheme, 
we can also show that any trajectory u{t) is separated from the singularities ±1 if, in 







lim F{u) = oo, F{u) : 
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addition, 



fin) 



> 



C 



p> 1 



u 




(see Remark 4.9; see also [8j for a similar condition for the Caginalp system). In that 
case, we have |M(t,x)| < 1 — 5 for some 5 > and, consequently, the problem becomes 
factually nonsingular and can be further investigated by using the techniques devised for 
the Cahn-Hilliard equation with regular potentials. Unfortunately, this last condition is 
not satisfied by the physically relevant logarithmic potentials and we indeed need the 
variational inequalities (and solutions) in order to deal with such a potential. 

The next, natural, step is to study the asymptotic behavior of the system. In particular, 
we are interested here in the study of finite-dimensional global attractors. We recall that 
the global attractor is the smallest compact set of the phase space which is fully invariant 
by the flow and attracts the bounded sets of initial data as time goes to infinity; it 
thus appears as a suitable object in view of the study of the longtime behavior of the 
problem. Furthermore, when the global attractor has finite dimension (in the sense of 
covering dimensions such as the fractal and the HausdorfF dimensions) , then, even though 
the initial phase space is infinite-dimensional, the dynamics of the system is, in some 
proper sense, finite-dimensional and can be described by a finite number of parameters. 
We refer the reader to, e.g., [2], [IH] and the references therein for extensive reviews 
and discussions on this subject. One powerful method, in order to prove the existence of 
the finite-dimensional global attractor, is to prove the existence of a so-called exponential 
attractor (in particular, this approach does not necessitate, contrary to the usual one, 
based on the Lyapunov exponents, the differentiability of the underlying semigroup). 
An exponential attractor is a compact and semiinvariant set which contains the global 
attractor, has finite fractal dimension and attracts all bounded sets of initial data at an 
exponential rate. We refer the reader to, e.g., [11], [13], [Hj and [36] for more details and 
discussions on exponential attractors. 

We thus prove the existence of global and exponential attractors for our problem. 
We emphasize that this result is obtained under general assumptions (without any sign 
assumption or any assumption of the form fll.7p ) and is thus valid for the variational 
solutions (which may not be solutions in the usual sense). In particular, such solutions 
may reach the singularities ±1 on sets of positive measure on the boundary x F or 
even on the whole boundary M"*" x F. This fact does not allow us to use the techniques 
devised in [3l| to establish the existence of finite-dimensional attractors for the singular 
Cahn-Hilliard system with usual boundary conditions (these techniques are strongly based 
on the fact that uit) is separated from the singularities for almost alH > 0, which is not 
true in our case in general). Instead, we prove the finite-dimensionality of the global 
attractor by using a proper modification of the techniques developed in [15j for porous 
media equations. 

This paper is organized as follows. In Section 2, we define proper (regular) approxima- 
tions of the singular potential and derive uniform (with respect to these approximations) 
a priori estimates which allow us, in Section 3, to formulate the variational inequality as- 
sociated with the singular Cahn-Hilliard system with dynamic boundary conditions and 
verify the existence and uniqueness of a solution for this inequality. We also study the 
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further regularity of the solutions. Then, in Section 4, we give sufficient conditions which 
ensure that the solutions are separated from the singularities of / and, thus, satisfy the 
equations in the usual (distribution) sense. Section 5 is devoted to the asymptotic be- 
havior of the system. Finally, we give, in Appendix 1, several auxiliary results. We also 
construct a simple example which shows that the solutions may not satisfy the dynamic 
boundary conditions in the usual sense for logarithmic potentials. 

2. Approximations and uniform a priori estimates 
We consider the following equations (for simplicity, we set all constants equal to 1): 



(2.1) 



dtU = A^fi, dnfi\^ = 0, 
fi = -A^u + f{u) + hi, m|^^q = Mo, 



in a bounded smooth domain Q of R^, endowed with dynamic boundary conditions on 

r := dn, 



(2.4) 



(2.2) dt-^ - Ari/J + giip) + dnu = h2, := u\ 

Here, u and fi are unknown functions, A^. and A-p are the Laplace and Laplace-Beltrami 
operators on Q and F, respectively, / and g are known nonlinearities, hi G L'^{Q) and 
/i2 £ -^^(r) are given external forces and 9„ stands for the normal derivative, n being the 
unit outer normal to F. 

We assume that the nonlinearity / has the form 

(2.3) f{z) := f{z) - Xz, 

where A G M is a given constant and the singular function / satisfies 

fl. /gC2((-1,1)), 

2. /(0)=0, lim„^±i /(n) = ±oo, 

3. /'(«)> 0, lim„_±i /'(m) = +00, 
^4. sgnw ■ f"{u) > 0. 

Since the function / is defined on the interval (—1, 1) only and has singularities at ±1, 
we a priori assume that 

(2.5) |M(t,x)| < 1 almost everywhere in M"*" x Q. 

We finally assume that the second nonlinearity g is regular on the segment [—1, 1], 

(2.6) geC\[-l,l]). 

Then, we can assume, without loss of generality, that g is smoothly extended to the whole 
line, g G C^(M), and g{z) = z + go{z) with ||5'o||c2(k) < C for some positive constant C. 

In order to solve the singular problem fl2.ll) . we approximate the nonlinearity / by the 
following family of smooth functions: 

(f{u), \u\<l-l/N, 

(2.7) fM := I /(I - l/N) + /'(I - 1/N){u - 1 + 1/N), u>l- 1/N, 

I /(-I + l/N) + /'(-I + l/N){u + 1 - l/N), u<-l + l/N, 
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and we set /Ar(M) := fNiu) — \u. We then consider the approximate problems 

\lJi = -^xU + fN{u) + hi, m|^^q = mo, 

endowed with the same dynamic boundary conditions (12. 2p . 

The main aim of the present section is to derive several uniform (with respect to — > 
oo) a priori estimates for the solutions (m, /i) = (mat, //at) of problems (12.81) . (12. 2p which will 
allow us (in the next section) to pass to the limit N —>■ oo and establish the existence of 
a solution for the singular problem (the existence, uniqueness and regularity of solutions 
for the regular case, such as in the approximate problems (12.81) . (12.21) . are now well- 
understood and will not be considered in the present paper, see |24j, [25j, |35], |12] and 
[is] for detailed expositions and related problems). 

As usual, it is convenient to rewrite problem (12.81) in an equivalent form by using the 
inverse Laplacian A := (— A^.)~^ (endowed with Neumann boundary conditions). To 
be more precise, since the first eigenvalue of the Laplacian with Neumann boundary 
conditions vanishes, we assume that the operator A is defined on the functions with zero 
mean value only and maps them onto the functions with zero mean value as well. Then, 
applying this operator to both sides of (12.81) . we have 

(2.9) Adtu := {~A,)-^dtU = A,u - fN{u) - h + (fi) , 

where (v) stands for the mean value of the function v over Q. Furthermore, taking into 
account (12. 2p . we see that 

(2.10) ifi) = - {A^u) + {hiu)) + ih) = dt + {9{u))^ - (/i2)r " (hiu)) + {h) , 

where (f)p := J^ f-^v{x)dS. We also mention that problem (12. 8p possesses the mass 
conservation law 

(2.11) (n(t)) = («(0)) = c 

and, thus, (dtu) = and the left-hand side of (12. 9p is well-defined. Finally, keeping in 
mind the singular limit —>■ oo, we only consider the initial data mq for which c G (—1, 1). 
We start with the usual energy equality. 

Lemma 2.1. Let the above assumptions hold and let u be a sufficiently regular solution 
of (12.90 . Then, the following identity holds: 

(2.12) l{^\\VMt)\\hin) + IW^MmW) + {FNHt)), Ih + {huu{t))n+ 

+ (G(w(t)), l)r - ih,,uit))r) + Wdtuml-.^^^ + \\dMt)\\hir) = 0, 

where Fjy^z) := fiy{s)ds, G{z) := f^g{s)ds, {■,-)n and (■, ■)r stand for the inner 
products in L^(fi) and L'^{T), respectively, and := {Az,z)n. 

Indeed, multiplying (12. 9p by dtU, integrating over Q and by parts and taking into 
account (12. 2p . together with the identity {dtu) = 0, we deduce (12.120 . 
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Corollary 2.2. Let the above assumptions hold and let, in addition, N be large enough. 
Then, any (sufficiently regular) solution u of problem f l2.9p satisfies: 

(2.13) + ||w(t)||^,(r) + {FN{u{t)), 

+ / mn{s)\\l-.^^) + \\dtu{s)\\l.^^^)ds< 
Jo 

<c{\\um 

where Ff^[z) := JJ^ fN^s) ds and the constant C is independent of t and m(0). 

Indeed, owing to our assumptions on / and the explicit form of the approximations /jv, 
see ( 12 ■71) . we can easily show that 

(2.14) 2Fn{z) + C> Fn{z) > ^Fr,{z) - C 

if N > Nq{X) is large enough, where the constant C only depends on A. Integrating now 
(I2.12P with respect to t and using fl2.14p . the fact that go{u) is globally bounded and 
obvious estimates, we end up with fl2.13p . 

As a next step, we obtain the dissipative analogue of estimate fl2.13p . 

Lemma 2.3. Let the assumptions of Lemma \2. 1\ hold, u be a sufficiently regular solution 
of (12. 9p and N be large enough (depending on A and c = (uq)). Then, the following 
estimate holds: 

(2.15) ||M(t)||^,(^) + ||n(t)||i,(p) + (Fjv(n(t)), 1)^+ 

+ I ^ (115*^(^)11^-1(1.) + \\dMs)\\l,^^) + ||/^(«(s))|Ui(o)) ds < 

< ciWummm + hmimiT) + (FMO)), i)n)e-"* + C(l + \\h,\\l,^^^ + Wh^Wl.^^^), 

where the positive constants C and a are independent of N and u, but can depend on the 
value c in the mass conservation (12.111) . In addition, the following smoothing property 
holds: 

(2.16) ||w(t)||^i(f,) + |k(t)||^i(r) + {F^{u{t)),l)n < 

< Ct-\\\u{Q) - c||^_,(^) + Wumhiv) + WhiWhi^n) + I|/^2||i2(r) + 1), t G (0, 1], 
where the constant C is independent of N. 

Proof. We have, owing to assumptions (12. 4p on the nonlinearity / and the fact that 

ce (-1,1), 

(2.17) fN{z).{z-c)>afN{z)z-C>a/2\fN{z)\~Cu z eR, 

where N is large enough and the positive constants a and Ci depend on c and A, but are 
independent of N (see |35]). Multiplying now equation (12.91) by u{t) := u(t) —c and using 
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the above inequality, we find 

(2.18) + + 

+ MtWHiin) + Mt)\\miT)) < C(l + WhWh^^) + I|/i2|li2(r)), 

for some positive constants a and C. Applying the Gronwall inequality to this relation, 
we obtain 

(2.19) \\u{t)\\l-Hn) + Mt)\\hir) + 

+ (ll^(s)llii(n) + ll^(s)llHi(r) + (/^(«(s)>w(s))n)c^s < 

< C{\\m\\l-^in) + Il^(0)||i2(r))e-"* + C(l + + ||/i2||i2(r)) 

for some positive constants C and a. In order to finish the proof of the lemma, there only 
remains to note that, owing to the monotonicity of the function /^r, 

(2.20) Fn{z) < fN{z).z, zeR. 

Then, the smoothing property fl2.16p follows in a standard way from 02.131) . 02.191) and 
O2.20p and the dissipative estimate 02.15P is an immediate consequence of the dissipative 
estimate 02.190 (in a weaker norm) and the smoothing property 02.16p . together with 
02.13p . This finishes the proof of Lemma [2. 3[ □ 

We are now ready to obtain additional regularity on dtu(t). To this end, we differentiate 
equation 02. 9p with respect to t and set 9(t) := dtu{t). Then, this function solves 

(2.21) {-A^)-^^te = AJ~fM^ + {^t^l), el^^ = 9,, 

where := -A^(Aj,Mo - /Af(Mo) - hi), and 

dtO - Ar^ + dnO + g'{u)e = 0, on T. 

Lemma 2.4. Let the assumptions of Lemma \2.1\ hold. Then, the following estimate is 
valid for the derivative 6{t) := dtu{t): 

(2.22) \m\\l-Hn) + ll^Wlli^(r) + l'^\ms)\\l.^n) + ds < 

< C(h(0)||i.(^) + ||n(0)||^.(r) + imfH-Hn) + 11^(0) ||i.(r))e-°*+ 

+ C{1 + ||/il||L2(f1) + ll^2||L2(r))' 

where the positive constants C and a can depend on the total mass c, hut are independent 
of N. In addition, the following smoothing property holds: 

(2.23) ||^(t)||?,-i(n) + ll^(t)lli2(r)< 

< Ct-'iMO) - c||^-,(^) + ||n(0)||i,(r) + + IMl^^r) + 1), (0, 1], 

where the constant C is independent of N. 
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Proof. We multiply equation fl2.2ip by 6{t), integrate over Q and use the fact that 
fN^'^) — Then, using also the boundary conditions and the fact that g'{u) is uniformly 
bounded, we find 

(2.24) j^iimWH-Hn) + ll^Wlli.(r)) + "(II^W llii(n) + imWrnir)) < 

<Cmumhin) + \\dMt)\\hir)), 

for some positive constants a and C which are independent of A^. Interpolating between 
and and applying the Gronwall inequality to this relation, we obtain the desired 
estimate (12.221) . Combining this estimate with (12.151) and (12.161) and arguing in a standard 
way, we end up with (I2.23P and finish the proof of the lemma. □ 

The next lemma gives i/^-estimates on the solutions for every fixed time t > 0. 

Lemma 2.5. Let the above assumptions hold. Then, for every fixed t > 0, the following 
estimate holds: 

(2.25) ||M(t)||^,(^) + l|M(t)||^i(p) + WfNiumiLHu) < 

< C(l + \\dMt)\\l-iin) + WdtumhiD + + Il^2||i2(r)), 
where the constant C depends on c, hut is independent oft and N . 

Indeed, multiplying equation (12. 9p by u{t) := u{t) — c and arguing as in the derivation 
of (I2.18P (but now without integrating with respect to t), we deduce the desired estimate 
(I2.25p . Here, we have used the inequality (12.170 again. 

Furthermore, using (12.250 and expression (I2.10p for the mean value of fi, we have 

(2.26) I (/i(t)) I < C(l + \\dMt)\\l-^(n) + WdtumhiD + IMl^n) + IMl^r))- 
We finally rewrite equation (12.91) in the form of a nonlinear elliptic problem, 

( A,u(t) - fN{u{t)) - u{t) = hit) := hi- 

(2.27) } -u{t)-Xu{t) + {-A,y'dtu{t)-{fi{t)), inQ, 
(Aru{t) - u{t) - dnu{t) = h2{t) := -hi + QoHt)) + dtu{t), on T, 

for every fixed t and note that the estimates derived above yield the following control of 
the right-hand sides in (I2.27P : 

(2.28) ||/ii(t)|U2(n) + ||/^2(t)||L2(r) < 

< C(l + ||5i«(t)||l,-i(f,) + ||atw(t)||i.(r) + + WWl.^r)) 

for some positive constant C which is independent of A^. 

Therefore, additional smoothness on the solution u := un can be obtained by a proper 
elliptic regularity theorem (see [35]). Unfortunately, in contrast to the case of regular 
potentials, this problem does not satisfy the maximal regularity estimate in for singular 
potentials /, see Appendix 1. Nevertheless, the partial regularity formulated below is 
crucial for what follows. 



10 



A. MIRANVILLE AND S. ZELIK 



Lemma 2.6. Let the above assumptions hold and set fi^ := {x G Q, d{x, T) > e}. Denote 
by n = n{x) some smooth extension of the unit normal vector field at the boundary inside 
the domain Q. Let also DrU := V^m — {dnu)n be the tangential part of the gradient VxU- 
Then, for every e > 0, the following estimate holds: 

(2.29) + \\V,DMt)\\LHn) + ||^2(o,) + IhiWb^cr) < 

<CMHt)\\LHn) + \\h2{t)\M 
for some positive constants a and which are independent of N. 

The proof of this estimate is based on some variant of the nonhnear locahzation tech- 
nique and is given in Appendix 1 (see Theorem 16. II) . 

We summarize the a priori estimates obtained so far in the following theorem which is 
the main result of this section. 

Theorem 2.7. Let the above assumptions hold and let u be a sufficiently regular solution 
of problem (12.91) with a sufficiently large N (depending on the constant A and the total 
mass c G (—1,1)^. Then, the following estimate is valid for every e > 0: 

(2.30) ||M(t)||^.(f,) + |k(t)||^.(r) + + \Ht)\\m(n) + 

+ \\dMt)\\H-Hn) + \\dMt)\\hir) + 

+ WV^DMrnlnn) + llfNiumWin) + J {\\dMs)\\min) + \\dtu{s)\\l,^^^) ds < 

< C{1 + Wurn'mm + hmimir) + WdMOWn-Hn) + 11^*^(0) |li.(r))'e-^*+ 

+ C(l+||/ii||i.(^) + ||/.2||i.(r))', 

where the positive constants a, [3 and C (which can depend on e) are independent of 
N ^ oo. In addition, the following smoothing property holds: 

(2.31) \\dtu{t)\\H^.^^) + \\dtu{t)\y^r) < Ct-W\u{{)) - c\\H-^n) + ||^i(0) |U2(r) + 

+ ||^i||L2(n) + ||^2||L2(r) + 1), i e (0, 1], 
where the constant C is also uniform with respect to N oo. 

Remark 2.8. We thus have a uniform if ^-estimate on the solution u inside the domain 
and an L^-estimate on the gradient of the tangential derivatives VxDtU. In contrast to 
this, we do not have a uniform control of the second normal derivative d^u close to the 
boundary. Nevertheless, since the L^-norm of the nonlinearity /at is controlled, (12.271) . 
together with the control of the tangential derivatives, allow us to estimate the L^-norm 
of d'^u. We thus have the control 

(2.32) \\u{t)\\w2,i(Q) < C{1 + \\hi\\L^n) + WhWL^r))- 

This, in turn, gives a control of the L^-norm of the normal derivative dnU at the boundary 
(owing to a proper trace theorem), 

(2.33) \\dnu{t)\\Li(^r)<C\\u{t)\\ 
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As we will see in the next section, estimates f l2.32p and fl2.33p remain true for the limit 
(as —>■ oo) solution u of the singular problem as well and, consequently, the trace of 
dnu{t) at the boundary is well-defined. However, owing to the nonreflexivity of L^-spaces, 
this trace may not coincide with the limit of dnUM{t)\^ computed on the boundary by 
using the dynamic boundary condition (12. 2p . Therefore, the boundary condition (I2.2p 
may be violated for the limit singular solution. As we will see below, this indeed happens, 
even in the ID case with smooth data. We overcome this difficulty by using monotonicity 
arguments and a proper variational formulation of problem (12. 9p . see Section [31 

We conclude this section by establishing the standard uniform Lipschitz continuity of 
the solution u of problem (12.90 with respect to the initial data. 

Proposition 2.9. Let the above assumptions hold and let Ui{t) and U2{t) be two (sujfi- 
ciently regular) solutions of problem (12. 9p such that 

{uim = {u,m = c. 

Then, the following estimate holds: 

(2.34) ||Mi(t) - U2{t) llH-i(n) + \\uiit) - u^it) |U2(r) < 

< C{\\ui{Q)-U2m\H-Hn) + ||mi(0) -M2(0)|U2(r))e^*, 
where the constants C and K are independent of t, N, Ui and U2- 
Proof. Let v{t) = Ui(t) — U2{t). Then, this function solves 

dtv - Arv + dnV + [g{ui) - g{u2)] = 0, on T. 

Taking the scalar product of the first equation with v(t), integrating by parts and using 
the facts that {v(t)) = 0, f'j^{z) > —X and the nonlinearity g' is globally bounded, we 
obtain 

dt 



(2.35) 



(2.36) |(||t;(t)||^.(^) + ||t;(t)||i.(r)) + 



for some positive constants a and C which are independent of A^. Interpolating between 
and in order to estimate the L^-norm of f in f2 and applying the Gronwall 
inequality, we find (12.340 and finish the proof of the proposition. □ 

3. The singular problem: variational formulation and well-posedness 

The aim of this section is to pass to the limit A^ — oo in (12. 9p and prove the existence 
and uniqueness of solutions of the limit singular problem (12.10 . As already mentioned, 
this limit solution is not necessarily a usual distribution solution of the equations and we 
need to define it in a proper way. To this end, we first fix a constant L > such that 

(3-1) l|VxM||i2(f^) - \\\u\\l2^^^ + L\\u\\l-i^^) > l/2||M||^i(f2) 

for all u G H^{Q) with (u) = and introduce the quadratic form 

(3.2) B{u, v) := (V,M, V,v)n - A(m, v)n + L((-A,)~'m, v)n + (Vrw, Vrt')r, 
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where u,v E H^{Q) fl -f^^(r) := {u, u G H^{VL)^ u\y € H^{T)} (and a similar definition 
holds for similar spaces) and u := u — (u) , v = v — (v) . Then, obviously, 

(3.3) B{u,u) = B{u,u) >0 

for all u G H\n)nH\T). 

The limit problem (I2.9p . corresponding to = oo, reads, formally, 

' {—Ax)~^dtU = AxU — f{u) + Am + (/x) — hi, in f2, 
^g^^ ^ H:= -AxU + f{u) - \u + hi, m|p = ^, 

(9t?/' - ArV" + fi'(^) + <9„M = /i2, on r. 

We test the first equation (again formally) with the function u — v, where v = v{t,x) is 
smooth and satisfies 

(n(t) - v{t)) ^ 0. 
Then, after an integration by parts, we have 

{Adtu, u-v)n + {dtu, u - v)y + B{u, u - v) + {f{u), u - v)n = 

= L{Au, u-v)n- {g{u),u - v)r - (h, u-v)n + (/i2, u - v)r. 
Finally, since B is positive and / is monotone, we have 

B{u, u-v)> B{v, u-v), {f{u), u-v)n> {f{v),u - v)ii, 

which yields 

{AdtU, u-v)n + {dtU, u - v)r + B{v, u - v) + {f{v),u - v)n < 

< L{Au, u-v)n- {g{u), u - v)r - (h, u-v)n + (/la, u - v)r. 

We recall that this inequality holds (again formally) for any properly chosen test function 
V such that (f (t)) = c. We are now ready to define a variational solution of the limit 
problem (13.41) . 

Definition 3.1. Let 

(3.5) (no,^o)G<f :={(^/,^)GL°°(fi)nL°°(r), \\u\\ ^n) < I , ||^||Loo(r) < 1}. 

A pair of functions {u, ip), u = u{t, x), x & Q, ip = tp{t, x), x G F, is a variational solution 
of problem (13. 4p if 

(3.6) ^(^)|r ~ ^(^) almost all t > 0, u{0) = uq, "ipiO) = V'o, 

1) - 1 < u{t,x) < 1 for almost all (t,x) G x 1], 

(3.7) 2) {u,iP) e C{[0,oo),H-\n) X L\T)) n L^{[0,T],H\n) X H\T)), VT > 0, 
3) fiu) G L\[0,T] X n), idtu,dt^) G L^i[T,T], H-\n) x L''{T)), VT > r > 0, 

{u{t)) = ('u(O)) and the variational inequality 

(3.8) {Adtu{t),u{t) - w)^ + {dtu{t),u{t) - w)v + B{v, u{t) - w) + (/(w), u{t) - w)n < 
< L{Au{t),u{t) - w)^ - {g{u{t)),u{t) - w)r - ihi,u{t) - w)n + (/i2, u{t) - w)r 
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is satisfied for almost every t > and every test function w = w{x) such that 

w e H\n) n H\T), f{w) e L\n) 

and (w) = (u(0)). Note that the relation u(t)|p = 'ip{t) is assumed to hold only for t > 0. 
At the initial time t = 0, no relation between uq and ipo is assumed. However, for t > 
the function can be found if u is known. Therefore, we can indeed write the variational 
inequality (13.81) in terms of the function u only. 

Before studying the existence and uniqueness of variational solutions, it is convenient 
to rewrite the variational inequality in terms of test functions v = v{t,x) depending on 
t and X. More precisely, let the test function v satisfy the regularity assumptions (13.71) 
and {v{t)) = (uq) = c (we will call this class of functions admissible test functions below) . 
Then, we can write inequality (13.81) with w = v{t) for almost all t > 0. Moreover, due 
to the regularity assumptions (13.71) on u and v, we see that all terms obtained are in 
with respect to t. Thus, we can integrate this inequality with respect to t, which gives 

(3.9) / [{Adtu,u-v)Q + {dtU,u-v)r]dT+ 

J s 

+ ^ [B{v, u-v) + {f{v),u - vU dr < 
< / [L{u,A{u-v))n-{g{u),u-v)r-{hi,u-v)n + {h2,u-v)r]dT 

J s 

for alH > s > 0. 

The next theorem gives the uniqueness of such variational solutions. 

Theorem 3.2. Let the nonlinearities f and g and the external forces hi and h2 satisfy 
the assumptions of Section Then, the variational solution of problem (13.41) (in the 
sense of Definition \3. 1\ ) is unique and is independent of the choice of L satisfying (13.11) . 
Furthermore, for every two variational solutions ui and U2 such that (mi(0)) = (^2(0)), 
the following estimate holds: 

(3.10) \\ui{t) - U2{t)\\H-^(^n) + Ui{t)-Mt)\\L\T) < 

< Ce^*(||wi(0) - u^miH-Hn) + ||^i(0) - ^^2(0) |U^(r)), 
where the constants C and K are independent of t, Ui and U2- 

Proof. We first need to deduce one more variational inequality for a solution u. Let w be 
a test function satisfying the assumptions of Definition 13.11 and set 

Va := (1 — Oi)u + aw, a G [0, 1]. 

Then, owing to assumption (I2.4p (4). the function \f{u) \ is convex and, therefore, 

\f{v^)\<\f{u)\ + \f{w)\. 

Consequently, Va is an admissible test function for every a G [0, 1]. Inserting f = f q in the 
variational inequality (13.91) . dividing it by a and using the fact that u G AC{\s, t], if~^(f2)n 
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L^(r2)) (here, AC stands for absolutely continuous), we see that 
(3.11) / [{Adtu,u - w)n + idtu,u - w)r]dT+ 

J s 

+ / [B{vo,,u - w) + {f{va),u - w)n]dT < 

J s 

< / [L{u, A{u - w))n - igiu),u - w)r - {hi,u - w)n + {h2,u - w)r]dT. 



Passing to the limit a — > in (13.111) and using the Lebesgue dominated convergence theo- 
rem for the nonlinear term, we end up with the desired additional vartiational inequality, 
namely, 

(3.12) / [{AdtU,u-w)n + {dtu,u-w)T]dT+ 

J s 

+ / [B{u,u - w) + {f{u),u - w)n]dT < 



< / [L{u,A{u-w))n- {g{u),u-w)r- ihi,u-w)n + {h2,u-w)r]dT, 



where w = w{t,x) is an arbitrary admissible test function. 

We are now ready to prove the uniqueness. Let Ui and U2 be two variational solutions 
of problem ( 13. 4p . We consider the variational inequality ( 13. 9p with u = ui and v = U2, 
together with the additional variational inequality f l3.12p with u = U2 and w = ui (this 
makes sense, since ui and U2 are admissible test functions), and sum the two resulting 
inequalities. Then, the terms containing 5, /, hi and /i2 vanish and, using, in addition, 
the fact that Ui G AC{[s,t\,H''^{Vt) fl L'^{T)), i = 1,2, we end up with the following 
inequality: 

(3.13) ^(||ni(t) - n2(t)||^-i(f2)nL2(r) " WM^) - U2{sWH-i(n)nLHr)) < 



< 



[L||ni(r) - M2(r) - (giM^)) - ^(«2(r)), Mi(r) - n2(r))r] dr. 



where ||'ii||jf--i(n)nL2(r) = + ll'"llL2(r)- Using now the fact that g G C^([— 1,1]) 

and applying the Gronwall inequality to f l3.13p . we see that 

\\ui{t) -U2(t)||i-i(n)nL2(r) < C'e^^*~"^lki(s) - «2(s)||H-i(n)nL2(r) 

for some positive constants C and K which are independent of t > s > and ut, i = 1,2. 
Passing to the limit s ^ in this estimate and using the continuity (I3.7p (2) of Ui and U2, 
we deduce the desired estimate fl3.10p which, in particular, gives the uniqueness. 

Thus, we only need to prove that the above definition of a solution is independent of 
the choice of L. To this end, we assume that ui is a variational solution for L = Li and 
U2 is a variational solution for L = L2. Let also Ui{0) = ^2(0). Using then the obvious 
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relation 

Bl^{v,ui -v) -Li{ui,A{ui -v))n = 

= Bl^{v,ui -v) - L2{ui,A{ui - v))n - (Li - L2)\\ui - 

and arguing exactly as in the proof of (13.101) . we have 

(3.14) ^{\\ui{t) - M2(t)||i-i(n)nL2(r) - WMs) - «2(s) ||^-i(f7)nL2(r)) < 



< 



^ [Li\\Ui{t) - M2(r)||^-i(f^) - {g{Ui{T)) - g{u2{T)),Ui{T) -U2{T))r]dT, 



which coincides with (13.131) and, therefore, also leads to estimate (I3.10p . Thus, Ui = U2 
and Theorem 13.21 is proved. □ 

We are now able to prove the existence of a variational solution u of problem (13. 4p by 
passing to the limit — > cxd in equations (12.91) . 

Theorem 3.3. Let the assumptions of the previous theorem hold. Then, for every pair 
{uo,ipo) G problem (13.41) possesses a unique variational solution {u,ip) in the sense 
of Definition \3.1\ Furthermore, this solution regularizes as t > and all the uniform 
estimates obtained in Section hold for the solutions of the the limit singular equation 
(13. 4p . In particular, the following estimate is valid for every e > 0: 

(3.15) |k(t)||2.„(^) + ||^i(t)||^,(r) + 
+ \Ht)\\l2^^^) + \\u{t)\\ 

+ Wdtuml-Hn) + \\dtumhir) + 
i-t+i 

+ \\V.D,u{t)\\hin) + \\f{umW(n) + mu{s)\\l.^n) + \\dtu{s)\\l.^r)) ds < 

,t' + l 



for some positive constants a and C which are independent of t and u (we recall that 
fl^ := {x G fl, d{x,T) > e}), where DrU denotes the tangential part ofVxU (see Lemma 
2.6); in addition, all norms in the left-hand side of (I3.15p make sense for any variational 
solution u. 

Proof. Let mat be the solution of the approximate problem (12. 9p . Then, repeating the 
derivation of the variational inequality (13. 9p . we see that 

(3.16) / [{AdtUN.UN -v)^ + {dtUN,UN -v)Y\dT+ 

J s 

+ {B{v, un-v) + {fN{v),UN - v)n) dr < 
< / [L{uN,A{uN-v))n-{g{uN),UN-v)r-{hi,UN-v)n + {h2,UN-v)r]dT 
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for every admissible test function v and every t > s > (we recall that the solution ujv of 
the regularized problem fl2.9p is smooth and all the formal calculations performed in the 
derivation of (13. 8p can be easily justified in that case). 

Our aim is to pass to the limit — oo in (I3.16p . We start with the case where the 
initial datum uq is smooth and satisfies the additional conditions 

(3.17) \uo{x)\<l-6, 5>0, ^o:=%|r. 

Then, according to Theorem 12. 7[ the sequence un satisfies the uniform estimate (12.301) 
and, therefore, we can assume, without loss of generality, that un converges to some limit 
function u in the following sense: 

1) un-^u weakly-* in L~([o,r], (i^^(^])^i^2 (r))n/f2(fi,) ), 

2) dtUN dtu weakly-* in L°°([0, T], //-^(r]) n ^^(r)) 

(3.18) and weakly in L^{[Q,Tl H^{VL) n H^{T)), 

3) DluN ^ Dlu weakly-* in L°°([0, T], ^^(fi)), 

4) u strongly in C'^([0,T] x Q) for some 7 > 0. 

Indeed, the last strong convergence follows from the facts that mat is uniformly bounded 
in L°°{[0,T],C°{n)), a > 0, and ^iM^v is uniformly bounded in L^{[0,T], H-\n)) (owing 
to Theorem 12.71 and the assumption that the initial datum uq is smooth and is separated 
from the singularities ±1). 

These convergence results allow us to pass to the limit — cxo in (13.161) and prove 
that the limit function satisfies (13. 9p for any admissible function v. The only nontrivial 
term containing the nonlinearity /at can be treated by using the inequality 

\fN{v)\ < \f{v)l 

the fact that f{v) G L^([0,T] x Q) and the Lebesgue dominated convergence theorem. 

Thus, we only need to show that the function u thus constructed satisfies the regularity 
assumptions (13.71) . The only nontrivial statements that we need to prove are that (I3.7p (l) 
holds and f{u) G L^{[0, T] x fl) (the other ones are immediate consequences of (13.180 ). Let 
us check the first one. Since the L^-norm of /Ar(uAr) is uniformly bounded, we conclude 
from the expression of the function that 

(3.19) measiit, x) e[T,T + l]xn, \uM{t,x)\>l-^\ <^{h, M > N, 



where 

(3.20) !f{x 



C 



max{|/(l-a;)|,|/(x-l)|} 

for some constant C which is independent of T G M^, M > N and N eN. Thus, passing 
to the limit M, ^ 00 in (13.190 and using the fact that ip{x) — ^ as a; — 0, we conclude 
that 

meas{(t,x) G [r,T + 1] x Q, \u{t,x)\ = 1} = 

and (I3.7p (l) is verified. In order to prove that f{u) is integrable, there only remains to 
note that the already proved statement (13.70 (1). together with the convergence (13.180 (4). 
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imply the almost everywhere convergence /at (mat) — > f{u) and, therefore, owing to the 
Fatou lemma, 

(3.21) \\f{u)\\Li{[T,T+i]xn) < liminfAT^oo \\fN{uN)\\L^{[T,T+i]xn) < oo. 

Thus, u is indeed the desired variational solution of fl3.4p and estimate (13.151) immediately 
follows from fl3.18p and Theorem 12.71 (in order to deduce the L^-estimate on f{u), one 
needs to use, in addition, inequality (13.211) ). 

Finally, we are now able to remove assumption (13.171) . To this end, we approximate the 
initial datum (uq, ipo) G $ by a sequence (mq, V'o) smooth functions satisfying (13.171) (of 
course, with 6 = 6k which can tend to zero as A; ^ oo) in such a way that 

(3.22) \\uo-Uo\\L^{n) ll^^olr " ^oIIl^cd ^ 0, (wg) = (mq) . 

Let {uk{t),'ijjkit)) (where ipk = Uk\j.) be a sequence of variational solutions of problem 
(13. 4p satisfying {uk{0),ipk{0)) = (""o'V^o) (whose existence is proved above). Then, owing 
to the uniform Lipschitz continuity estimate (I3.10p and assumption (I3.22p . {uk,ipk) is a 
Cauchy sequence in C{[0,T], H~^{Q) x L^(r)) and, therefore, the limit function 

{u,ip) := lim {uk,ipk) 

exists and also belongs to C{[0,T], H-\T) x L'^{n)). The fact that M is a variational 
solution of (13. 4p . as well as estimate (I3.15p . can be verified, based on the uniform estimates 
derived in Section [2l exactly as was done above for smooth initial data. This finishes the 
proof of Theorem 13. 3[ □ 

Corollary 3.4. Under the assumptions of Theorem \3.3\ equation (13. 4p generates a solu- 
tion semigroup S{t) in the phase space $, 

(3.23) Sit){uo,iPo):={u{t),m), S{t) : <^ ^ <t>, t>0, 

where {u{t),ip{t)) is the unique variational solution of problem (13. 4p with initial datum 
{uq,iPq). Furthermore, this semigroup is globally Lipschitz continuous, 

(3.24) \\S{t)iuliPl)-Sit)iuliPl)\\H-^n)><LHr)<Ce''%ul-4^ 

in the metric of the space := H^^{^1) x L^(fi). 

We now start investigating the analytic structure of a solution {u{t),ip{t)) of problem 
(13. 4p (this will be continued in Section H]). 

Proposition 3.5. Let {u{t),ilj{t)) be a variational solution of problem (13. 4p constructed 
in Theorem \3.3l Then, il>{t) = u{t)\^ for t > and, for every ip G C^((0,T) x Vl) such 
that {^p{t)) = 0, there holds 

(3.25) / {i-A,)-'dtuit),ipit))ndt = 
Jk+ 

= I ((A,.n(t), <^(t))n - (/(u(t)), v{t))n + A(u(t), ^{t))n - (/^i, <^(t))n) dt. 
Jr+ 

Furthermore, 



(3.26) 



M e L°°([r,r],W^2'^(fi)), T>r>0, 
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and the trace of the normal derivative on the boundary, 

(3.27) [d^uU:=d^u\^eL'^i[T,T],L\r)), T > r > 0, 
exists. 

Proof. Since, according to Theorem 12.71 the approximating sequence u^- is uniformly 
bounded in L'^{[t,T], H'^{fl^)), for any e > 0, the sequence fNiu^) is also uniformly 
bounded in L'^{[t,T], L'^{Qi;)). This fact, together with the almost everywhere conver- 
gence established in the proof of Theorem 13. 3[ guarantee that /at (mat) f{u) weakly in 
L^([r, T] X Vie) for all e > and this, in turn, allows us to verify identity (13.251) by passing 
to the limit in the analogous identity for the approximate solutions ujsf. 

In order to check the remaining statements of the proposition, we first deduce from 
(1^ that 

(3.28) {-^^Y^dtu{t) = AMt) - fW)) + \u{t) -hi + c{t) 

for some function c G L°°([r, T]), T > r > 0. At this point, equality (13.281) is un- 
derstood as an equality in ^^^^([r, T] x Vl). However, owing to estimate (13.151) . f{u) G 
L^{\t,T], L^{VL)) and the term (— A^;)^"'^^*^ also belongs at least to this space. Thus, we 
see that 

A,«GL°°([r,T],Li(n)) 

and, consequently, since VxDt-u is controlled by (13.151) . we can finally conclude that (13.261) 
holds, which gives the existence of the trace (I3.27P and finishes the proof of Proposition 
ESI □ 

Note that, using the obvious fact that {{—/\x)^^dtu) = 0, we can find the explicit 
formula for the function c{t) in (13.281) . namely, 

(3.29) c(t) = (A,M(t) - f{uit)) + Xuit) -hi) = - (/i(t)) 

and, therefore, the first equation of (13. 4p is satisfied in a usual sense (say, as an equality 
in Lf^^{[T,T] X Q) or/and almost everywhere). 

We now investigate the third equation of (13. 4p (the equation on the boundary). Ac- 
cording to Theorem 12. 7[ we see that the approximating sequence {uNit),ipN{t)) satisfies 

\\dt'ipN{t)\\L°^{lT,T],L2{r)) + ||V'Af(^)||L2([r,r],H2(r)) < C 

and, therefore, using the fact that the approximate solutions satisfy the second equation 
of (13. 4p . we can assume, without loss of generality, that we have the convergence 

(3.30) [dnu]ext:= lira dnUNLeL^i[T,T],L\T)), T > r > 0, 

where the limit is understood as a weak-star limit in L°^([r, T], L^(r)). Then, obviously, 

(3.31) dtip - Arip + g{ip) + [dnu\ext = h2, on T, 

and, in order to verify that the variational solution (m, ^p) satisfies equations (13. 4p in the 
usual sense, there only remains to check that 

(3.32) [dnu]int = [dnu]ext for almost every {t,x) G x T. 
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However, as the example in Appendix 1 shows, this identity can be violated even in the 
simplest ID stationary case. In the next section, we formulate several sufficient conditions 
which ensure that (13.321) holds for every (variational) solution of (13. 4p . 

4. Additional regularity and separation from the singularities 

The main aim of this section is to study the analytic properties of the variational 
solutions u of problem (13. 4p . especially close to the singular points ±1. We start with the 
following result which gives an additional regularity on u{t,x) close to the points where 
\u{t,x)\ < 1. 

Proposition 4.1. Let the assumptions of Theorem \3.2\ hold and let u be a variational 
solution of problem (13.41) . Let also 6 > 0, T > be given and set 

(4.1) ns{T):={xen, \u{T,x)\ <l-6}. 

Then, u G W'^''^{Vts{T)) and the following estimate holds: 

(4-2) \\u\\HHn,{T)) < Qs,T, 

where the constant Qs^t only depends on T and 6, but is independent of the concrete choice 
of the solution u. 

Proof. Since the solution m(T, x) is Holder continuous with respect to x (see (I3.15P ). there 
exists a smooth nonnegative cut-off function 9{x) such that 

r 1) e{x) = 1, xe ns{T), 

(4.3) h) 9{x) = 0, xen\ns/2{T), 

[3) ||6'||c2(K3) < Ks,T, 

where Ks^t depends on the constants in (I3.15p . but is independent of the concrete choice 
of the solution u. 

Furthermore, let UAr(t, x) be a sequence of approximate solutions of problems (12. 9p which 
converges to the variational solution u{t, x) as N 00. Then, since this convergence holds 
in the space C'^{[t,T] x Q) for some 7 > 0, 

(4.4) \un{T,x)\ < 1 - 6/4, X e ns/2{T) 

if is large enough. Set now vn{x) := 9{x)un(T, x). Then, this function obviously solves 

the following elliptic boundary value problem (compare with (12.270 ): 

(4.5) 

^A^Vn -vn = h.M := ef^MT)) + efiiiT) + 2VJ.V,UNiT) + un{T)AJ, 

^ Vn\j. = Wn, 

ArWN -wn- dnVN = h2{uN) ■■= 

:= 9h2{T) + 2Vr^.VrUjv(T) + M^(T)Ar0 - UN{T)dJ, 

where the functions hi, i = 1,2, are the same as in (12.271) . In addition, owing to estimates 
(1^ . (^M, <^M, <SM and (gaD, we see that 

(4.6) \\hiiuN)\\L2{n) + \\h2iuN)\\mr) < Q&,t, 
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where the constant Qs^t is independent of and of the concrete choice of the solution 
u. Applying the iJ^-regularity theorem to the linear elliptic problem fl4.5p (see P5|) and 
recalling (14.31) . we deduce that 

(4-7) \\uN{T)\\H2(^ns{T)) <Q5,T 

and, consequently, by passing to the limit — oo, we see that u{T) G H'^{Qs(T)) and 
(14. 2 p holds. This finishes the proof of the proposition. □ 

Remark 4.2. Applying the L^-regularity theorem to the elliptic boundary value problem 
(14. 5p . together with a proper interpolation inequality, we have 

\\''^{T)\\w'^'P(ns(T)) < Q5,p,T 

for any p < oo. However, it seems difficult to obtain further regularity results on u in 
^siT) by directly using equations (12.90 or (14. 5p . owing to the presence of the nonlocal 
term {fi{t)) (which is only with respect to t). Alternatively, one can use the standard 
interior estimates for the initial fourth-order problem (12. ip . Then, it is not difficult to see 
that the factual regularity of the solution u in fls(T) is only restricted by the regularity 
of the data /, g, hi, i = 1,2, and Q (and, if these data are of class C°°, the solution u is 
of class C°° in Qs{T) as well). 

Corollary 4.3. Let the assumptions of Theorem \3.2\ hold and let u be a variational solu- 
tion of problem (13. 4p . Assume, in addition, that 

\u(to,xo)\ < 1 

for some (to, Xq) G x T, with t^ > 0. Then, there exists a neighborhood {tQ — 6, tQ+6) x V 
of (to, a^o) M X r such that 

(4.8) [dnU]int{t, X) = [dnU]ext{t, x) , V(t, x) G (to - 6,to + 6) X V. 

In particular, if 

(4.9) |u(t,a;)|<l for almost all {t,x) eR'^ x T, 

then the equality [dnu]ext = [dnu]int holds almost everywhere in M"*" x F and, therefore, the 
variational solution u solves equations (13. 4p in the usual sense. 

Proof. Since the solution u is Holder continuous with respect to t and x, there exists 6 > 
such that the inequality 

\u{t,x)\ < 1 - (5 

holds for all (t, x) belonging to some neighborhood (to — 5, to + 5) x Vs of (to, xq) inRxfl. 
According to Proposition 14.11 the sequence un of approximate solutions (converging to 
the variational solution u) satisfies 

\\uN\\L°^{[to-5,to+5],H'2(Vs)) < C, 

where the constant C is independent of A^. Consequently, we can assume, without loss of 
generality, that un ^ u weakly-star in this space. Thus, 

weakly in L^([to — 6,to + 6] x V) (for a proper choice of the small neighborhood V of 
xq). This convergence, together with the definition (I3.30p of the function [dnu]ext, give 
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the desired equality (I4.8p . Thus, the first part of the statement is proved and the second 
one is an immediate consequence of the first one, which finishes the proof of Corollary 

K3\ □ 

Thus, in order to prove that any variational solution u of problem (13.41) satisfies the 
equations in the usual sense, it is sufficient to check (14.91) . The next corollary shows that 
this will be the case if the nonlinearity f{u) has sufficiently strong singularities at ±1. 

Corollary 4.4. Let the assumptions of Theorem \3.^ hold and let, in addition, the potential 
F{u) he such that 

(4.10) lim F{u) = oo. 

Then, for every variational solution u of problem (13.41) . 

F{u{t))eL\T) and \\F{um\LHT) < Ct 
for almost all t > T > and condition (14. 9 P holds. 

Proof. Let Mat be a sequence of approximate solutions converging to the variational so- 
lution u. Applying estimate (16. 4p in Appendix 1 to the elliptic problem (I2.27p . we infer 
that 

(4.11) WFNiuNmiLHD < Ct, t>T, 

where the constant Ct is independent of N. Using assumption (I4.10p and arguing as in 
the proof of Theorem \3.3\ we see that condition (14.91) indeed holds. Then, owing to the 
convergence ujv — > u in C^{[0,T] x fl), 7 > 0, we conclude that Fp^{uiy) F{u) almost 
everywhere in x F. The Fatou lemma finally yields that F{u{t)) G -^^^(r), which 
finishes the proof of the corollary. □ 

In particular, condition (14. 9 p is satisfied if the nonlinearity / is of the form 

(4.12) /(«) ~ 

with p > 1. Unfortunately, the assumptions of Corollary 14.41 are violated in the physically 
most relevant case of a logarithmic potential, 

1-1-7/ 

(4.13) /(^) = ln^. 

1 — u 

Furthermore, as explained in Appendix 1, in that case, the variational solution u may 
indeed not be a solution in the usual sense (even in the ID stationary case). However, the 
next proposition gives another type of sufficient condition (in terms of the nonlinearity 
g and the boundary external forces which guarantees the equality [dnu]ext = [dnu]int 
and holds for the logarithmic potential (14.130 . 

Proposition 4.5. Let the assumptions of Theorem hold and let, in addition, the 
following inequalities hold: 

(4.14) g{-l)+e<h2{x)<g{l)-e, xeT, 
for some e > 0. Then, condition (14. 9 p holds and 

(4.15) ||/(«)||Li([M+i]xr) < C,,T, t > T > 0, 
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where the constant Ci^^t is independent of the concrete choice of the variational solution 
u. In particular, every variational solution of fl3.4l) solves this system in the usual sense. 

Proof. As above, it is sufficient to derive the uniform (with respect to — oo) estimate 
( I4.15P for the approximate solution un oi (12.91) . In order to do so, we rewrite the system 
in the elhptic-parabohc form 

{^xUNif) - fN{uN{t)) - UN{t) = hi{t), 
dtipN - ArV'Af + dnUN + gii^N) = h2. 

Furthermore, since only the values of g on the segment [—1, 1] are important for the limit 
problem, we can assume, without loss of generality, that 

g{-u) + e < h2{x) < g{u) - e, ueR, |m| > 1, xeT. 

It follows from these inequalities and the continuity of g that 

(4.17) {g{z)-h2{x)).f^{z)>'-\fM{z)\ + Ce, zeR, xer, 

where the constant depends on e and g, but is independent of A^. 

We now multiply the ffist equation of (14.161) by /at (mat) and integrate with respect to 
X. Then, integrating by parts and using estimate (14.171) . we find 

(4.18) - / FN{uN{t)) dS + {f'j^{uN{t))V,UN{t),V^UN{t))n+ 

+ (/;(M^(t))VrU7v(t), Vrujv(t))r+ 
+ l/2||/^(n^(t))||i.(^) + 5/2||/^(n^(t))|Ui(r) < C(l + fhmhin))- 

Integrating this inequality with respect to time and using the facts that /j^ > and the 
L^-norm of hi(t) is controlled (see (I2.22p and (12.281) ). we obtain 

(4.19) ||/jv(^^^)||Li([M+i]xr) < '^{\\FN{uNm\LHr) + \\FN{uM{t+l))\\Li^r)) + Ce,T. 

There only remains to note that the right-hand side of (14.191) is controlled, owing to 
estimate (16.41) (exactly as in the proof of (14. lip ). Therefore, (I4.19P gives uniform bounds 
on the L^-norm of /at (mat) on the boundary. Passing to the limit N —>■ oo now gives the 
statement of the proposition (exactly as in the proof of Corollary 14.41) . □ 

Remark 4.6. As already mentioned, the variational solution u may not solve equations 
(13. 4p in the usual sense if conditions (I4.10p and (14.141) are violated (see Example 16.21 for 
details). Furthermore, arguing as in this example, it is not difficult to show that, for any 
singular nonlinearity / which does not satisfy (I4.10p . there exist "nonusual" variational 
solutions of problem (13. 4p if the external forces h2 are large enough. 

We conclude this section by estabhshing that every solution u of problem (13. 4p is sep- 
arated from the singularities ±1 if the nonlinearity / is singular enough. To this end, we 
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need to require at least condition (14.101) to be satisfied (see again Example 16.21) . Actu- 
ally, we will require slightly more, namely, that the nonlinearity / satisfies the following 
inequahties: 

(4-20) T^-kn < ^ < 



(1 -U2)P-1 - U - (1 -m2)A^ 

for some positive constants /tj, i = 1,2, and M and where p > 2 (recall that condition 
(I4.10p is violated if p < 2, so that the sufficient condition (I4.20p is close to the necessary 
one). In addition, we assume more regularity on the external forces hi and /i2, namely, 

(4.21) hi e L^{n), h2 e L°°(r). 

Theorem 4.7. Let the assumptions of Theorem \3.S\ hold and let, in addition, ( 14.20p and 

(I4.2ip be satisfied. Then, every variational solution u of problem (13. 4p is separated from 
the singularities ±1, namely, the following estimate holds: 

(4.22) \u{t,x)\<l-5T, t>T>0, xeQ, 
where the constant 6t depends on T, but is independent of u, t and x. 

Proof. We only give below the formal derivation of estimate (14.220 which can be justified 
as above by approximating the solution m by a sequence un of solutions of the regularized 
problem (12.90 . Our proof is based on the following lemma. 

Lemma 4.8. Let the assumptions of the theorem hold and let u(t) be a (variational) 
solution of problem (13.40 . Then, for every q > 0, f{u) G L'^([t,t + 1] x f2) for all t > 
and the following estimate holds: 

(4.23) \\fiu)\\Lmt,t+i]xn) < Ct,„ t>T>0, 
where the constant Cx^q is independent oft and u. 

Proof. We rewrite system (13. 4p in the form of a coupled elliptic-parabolic problem. 



(4.24) 



A^M - f{u) - u = hiit), u\^ = ip, 
dtip — A-pip + dnU = h2{t). 



Then, owing to the regularity estimate (I3.15P on the solution u and conditions (14.210 . we 
have 

(4.25) \\hi{t)\\mn) + Whrnio^iD < Ct, t>T. 

We introduce the function f{u) := Then, 

2u 

'1^2 



We multiply the first equation of (14.240 by u(p{u)^^^, where n > 1 is an arbitrary fixed 
exponent, and integrate with respect to x G fi. Then, integrating by parts and using the 
obvious transformations 

(v.n, vMi^(n)rv.«)n = ^(lv.w|^ i^'(w)iVHr"')n > c^niiv.(i^(w)r/2)iii.(^^ 
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and 



(owing to the first inequality in fl4.20p ). we have 



(4.26) ^||^(«)ir2.(r) + 2«:(|||^(«)r/l^.(^) + ||v^(n)||2tV)) < 
< C\\k{t)hs^n)Mu)\\ltL^^^^ 

where k > and the constant Ct depends on n, but is independent of t > T > and the 
concrete choice of the solution u. Let us estimate the right-hand side of (14.261) . In order 
to estimate the boundary term, we use the following trace inequality 

is+l ^ r^/llT/ll2 , IIT/Il2s 



V 



L(r) < + ll^ilL2''(n)) 



(which can be easily obtained by using a proper interpolation inequality and Sobolev's 
embedding theorem). Using this inequality with V = \^\"'^'^ and s = 1 +p/n, we find 

and, therefore, since p > 2, we can rewrite (14.261) without any boundary term in the 
right-hand side, 

(4.27) |||^(^)||^(r) + 2«:'(|b(t.)||23„(n) + \Mu)\\T4.^n)) + 

+ ^^i^^mtl^D < ^T(||v^(n)||^+U,)(^) + 1) 

for some positive constants k' and Ct which depend on n and T, but are independent of t 
and u (here, we have implicitly used the embedding C and replaced the exponent 
n-|-p/2byn-|-lin the boundary term). In order to estimate the right-hand side of this 
inequality, we use one more interpolation inequality, 

\Mu)\\Un) < lbMllt-(^)||(/.(^)||[l-/4) < C{Mu)\\ls.(^^ + Mu)\\lX^n))^ 
where s G [0, 1] and q are such that 

1 _ e 1-e 1 _ 9 1-e 

r 3n n + p' s n n + p 
Solving these equations for r = 3/2{n + 1) and n > p — 1, we have 

2Ti — T) 

^-("+l) 2n-p-(p-2) ^"+^ 

(since p > 2). Thus, 
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and yields, noting that p>2 and C L"", 

(4.28) |||^(«(t))||^(r) + ^'yHmitl.^n) + '^llv^(^W) < Ct- 

Multiplying this inequality by {t — T)"+^, we obtain 



(4.29) - Tr+'yHmur)] + ^'{t - Tr+'Mumit'.^n) < 

< -4it - T)Mum\LHr)r' + C{n + l)[(t - T)||v.(«(t))|U"(r)]"+ 

+ Crit - Tr+' < 

Integrating fl4.29p with respect to t G [T, T + 2] , we finally end up with 

Since n is arbitrary, this last inequality, together with the second inequality in (14.200 . 
finish the proof of the lemma. □ 

It is now not difficult to finish the proof of the theorem. To this end, we note that, 
owing to estimate (13.151) and Sobolev's embedding theorem, 

(4.30) ||M(t)|U2-i/3,3(r) < C\\u{t)\\H^r) <Ct, t> T. 

On the other hand, owing to the lemma and the first equation of (13. 4p . there holds 

\\^xU{t)\\Lg(^[t,t+l],L^{n)) < Ct- 

Thus, owing to the maximal regularity for the Laplacian in and Sobolev's embedding 
theorem, 

(4.31) \\VMt)\\LH[t,t+i]xn)<CT,g, t>T, 

for any q > 1. Furthermore, it follows from (13.151) . (I4.23P and (14.310 that 

IIV5('^^)l|L'-([t,t+l],VFi.'-(n)) < 

(4.32) < C(||v3(M)||L2r([i^(+i]xn) + ||v5'(m)||l2'-([M+1]xJ7))(1 + ||V:cM||L2r([i_i_^i]xf7)) < Cr,T, 

\\dtV{u)\\L'^-mt,t+i],L<^-^{n)) < C\\ip{u)\\Lrs{[t,t+i]xn)\\dtu\\L2([t^t+i]xH^in)) < C'e.r, 

where e > and r > are arbitrary and the constants Cr,T and Cs^t are independent of 
u and t > T. Fixing finally r ^ 1 and e ^ 1 in such a way that 

w^-'{[t,t + i] X n) n u'{[t,t + i],w^'''{n)) cc{[t,t + i] xn), 

we deduce from (14.321) that 

sup |- Y, ^1 = y{u)\\c{it,t+i]xn) <Ct, t> T, 



{s,x)(^[t,t+i]xn 



which gives (I4.22p and finishes the proof of Theorem 14. 7[ □ 
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Remark 4.9. It is not difficult to see that assumption (14.211) can be slightly relaxed, 
namely, 

hi G h2 G L"2(p)(r), ri(p) < 3, r2{p) < oo 

and ri{p) — > 3, r2(p) ^ oo as p ^ 2 (where p > 2 is the exponent in inequalities (I4.20p ). 
It is also worth noting that the proof of Lemma |48] does not involve the Laplace-Beltrami 
operator in the boundary conditions (and we have used it only in the derivation of estimate 
(I4.30p ). Consequently, arguing in a slightly more accurate way (e.g., by obtaining the L°°- 
estimate on f{u) directly from the estimates of Lemma l4.8l by a Moser iteration technique), 
one can extend the theorem to less regular boundary conditions, 

dtu + dnU + g{u) = /i2, X & V. 

Finally, the aforementioned Moser scheme also allows to remove the second inequality in 
(I4.20p and to obtain the separation from the singularities by only using the ffist inequality 
in (I4.20p for the function /. We will come back to these questions elsewhere. 

5. Long-time behavior: attractors and exponential attractors 

In this concluding section, we study the asymptotic behavior of the trajectories of the 
solution semigroup (13.231) acting on the phase space $, endowed with the metric of 
We ffist recall that problem (12.11) enjoys the mass conservation 

(5.1) Mt)) ^ KO)) := c. 

Therefore, it is natural to consider the restrictions of our semigroup to the hyperplanes 

(5.2) $,:={(«, V')G<I', {u) = c}, CG (-1,1), 5(t) : ^ 

The following proposition gives the existence of the global attractor Ac for this semigroup. 
We recall that, by definition, a set Ac C $c is the global attractor for the semigroup (15.21) 
if 

1) It is compact in $c- 

2) It is strictly invariant, i.e., S{t)Ac = Ac, t >0. 

3) It attracts $c as t — oo, i.e., for every neighborhood 0{Ac) of Ac in $c, there exists 
T = T{0) such that 

S{t)<^c^O{Ac), t>T. 
We refer the reader to, e.g., [2], [36] and [19] for details (we note that, in our situation, 
the phase space $c is, by definition, bounded and, therefore, we need not involve bounded 
sets in the definition of the global attractor). 

Proposition 5.1. Let the assumptions of Theorem \3.2\ hold. Then, for every c G (—1, 1), 
the semigroup S{t) associated with the variational solutions of problem (12. ip acting on 
the hyperplane (15. 2p (endowed with the metric of ) possesses the global attractor Ac- 
Furthermore, this attractor is bounded in the space C"(i7) x C"(r), a < 1/4, and is gen- 
erated by all complete trajectories of the semigroup S{t) (i.e., by all variational solutions 
{u{t),v{t)) which are defined for all t G Mj. 

Indeed, owing to estimate (I3.24p . the semigroup S(t) has a closed graph in <l>c- On the 
other hand, owing to estimate (13.151) . this semigroup possesses an absorbing set which is 
compact in $c (endowed with the metric of $"') and bounded in C"(f2) x C"(r). Thus, 
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the existence of Ac, together with all properties stated in the proposition, follow from a 
proper abstract attractor's existence theorem (see, e.g., [2], [SS] and [15]). 

Our next task is to prove the finite-dimensionality of the global attractor Ac constructed 
above and the existence of a so-called exponential attractor. We recall that, by definition, 
a set M.{c) C $c is an exponential attractor for the semigroup S{t) if 

1) It is compact in $c- 

2) It is semiinvariant, i.e., S{t)M{c) C M{c), t>0. 

3) It has finite fractal dimension in $c- 

4) It attracts $c exponentially fast as t — oo, i.e., 

dist$^(5(t)$e,A^(c)) < Ce-"", t > 0, 

for some positive constants 7 and C. Here and below, distv{X,Y) stands for the non- 
symmetric Hausdorff distance between sets in V. 

We also recall that the usual construction of an exponential attractor is based on the 
so-called squeezing (or smoothing) property for the difference of solutions (or their proper 
modifications, see tl3j, ^Mj and |36j for details). The main difficulty here lies in 
the singular nature of the equations at ±1. In particular, the difference between two 
singular solutions ui and U2 does not possess any regularization. Nevertheless, as we will 
see below, our nonlinearity is strictly monotone near the singularities ±1 and, far from 
these singularities, the problem still possesses the usual parabolic smoothing property. 
This fact, together with the Holder continuity of the solutions and some localization 
technique, allow to construct an exponential attractor by using a proper modification of 
the techniques developed in [13] and [TJ]. However, some additional difficulties arise here, 
due to the fact that the if~^-norm is not local. 

Theorem 5.2. Let the assumptions of Theorem \3.2\ hold. Then, the semigroup S{t) acting 
on the phase space <l>c (endowed with the metric of^) possesses an exponential attractor 
M{c) which IS hounded m C"(^]) x C°(r), a < 1/4. 

Proof. We first note that, owing to Theorem 13. 3[ there exists a compact (for the metric 
of $"') absorbing set 

(5.3) B, := 5(1)$, 
such that 

(5.4) S{t)Bc C Be 
and 

(5.5) ||M||c"([i,f+i]xn) + ||n(t)||H2(r) + \\dtu{t)\\H-^{Q.) + \\dtu{t)\\L'2(T) + 

+ ll/(^(^))l|Li(n) + \\dtU\\L^f^[t^t+l],m{n)) + \\dtU\\L-^{[t,t+l],H^(V)) < R 

for some fixed constant R which depends on c, but is independent of ('u(O), ^^(O)) G Be. 
In particular, for every point (m, ip) G Be, we have 

and, consequently, we generally write u instead of ip in the boundary norms. 

Thus, we only need to construct an exponential attractor M.{c) for the semigroup S(t) 
restricted to the semiinvariant absorbing set Be. As usual, to do so, we need to obtain 
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proper estimates on the difference of two solutions Ui{t) and U2{t) starting from tlie set 
Be- Furtliermore, we use tfie following natural norm on tlie pliase space B^: 



(5.6) \\ui - U2\\^w := \\ui - U2\\H-i{n) + ll^^i - «2||L2(r) 



2 

L2(r)- 



A crucial point in the proof is the global Lipschitz continuity in this norm (see Theorem 
2.2), 

(5.7) ||Ml(t) -M2(t)|||™ + 

t+1 

+ Iki(s) -M2(s)||^i(r))c?s < Ce'''\\u^{0) - U2{0)\\1^, 

where the positive constants C and K are independent of Mi(0),M2(0) G Be. 

We now consider an arbitrary small e-ball B{e, uq, $*") in the space Be (endowed with 
the metric of $*") and centered at uq, where < £ < £0 "C 1 (and the parameter Eq will 
be fixed below). Let also u^(t), t > 0, he the solution of problem (13.41) starting from uq- 

As in (14.11) . we introduce the sets 

(5.8) ^ls{uo) ■= {x e ^l, \uo{x)\ < 1 - 6} , Us{uq) := {x e ^, |uo(x)| > 1 - 5}, 

where 5 is a sufficiently small positive number. Then, since the function uo{x) is uniformly 
Holder continuous in Q, there holds 

(5.9) d{dns,{uo),dns^{uo)) > Cs,,52 > 0, 5i ^ ^2, 

where the constant C^^^^j depends on 6i, i = 1,2, but is independent of the concrete choice 
of Uo G Be. 

As a next step, we note that, owing to the uniform Holder continuity of the trajectory 
M°(t) (in space and time), there exists T = T{6) such that 

|n°(t)| <l-^, xeQsiuo), te[0,T], 
\u%t)\ > 1 - 2(5, xe n25{uo), t e [0,T], 

and, furthermore, owing again to the uniform Holder continuity, 

WMt) - Mt)\\cm < c\Mt)-u2mi>4Mt)-u2mc'4n) < Cts^ 

for all Mi(0),'U2(0) G B{e,Uo, $"'). We can thus fix Eq = eo{6) in such a way that 

^^^^^ \uit)\<l--, xensiuo), te[0,T], 

\u{t)\>l-46, xen25{uo), te[0,T], 

for all trajectories u{t) starting from the ball B{e,uo, $"") with e < sq. 
We also introduce the cut-off function 6 G C°°(R^, [0, 1]) such that 

(5.12) e{x) = o, xeTisiuo), e{x) = i, xen25{uo). 
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Such a function exists, owing to condition fl5.9p . Furthermore, it follows from this con- 
dition that this function can be chosen in such a way that it satisfies the additional 
assumption 

(5.13) ||^||c"=(iR^) - ^fc' 

where /c G N is arbitrary and the constant Ck depends on 5, but is independent of the 
choice of uq G Be, see |il5j for details. 

Finally the second estimate of fl5.1ip yields 

(5.14) f\u{t,x))>k{5), xeU^siuo), tG[0,T], 
for all trajectories u{t) starting from the ball B{e,UQ, $""), where 

A(5) :=min{/'(l-45),/'(-l + 45)}. 

Since f'{u) — > cxd as u — > ±1, then A{6) — >• cxd as 5 ^ and we can fix 5 > in such a 
way that A (5) is arbitrarily large. This will be essentially used in the next lemma which 
gives some kind of smoothing property for the difference of two solutions Ui and U2 and 
is crucial for our construction. 

Lemma 5.3. Let the above assumptions hold. Then, there exists 5 > such that the 
following estimate holds: 

(5.15) \\ui{T) - U2(T)|||„ < e-^^||Mi(0) - U2(0)|||.,+ 



T 



+ C \\e{u,{s)-U2{smUn)ds, 
Jo 

where the positive constants 5, P and C are independent ofui,U2 G B{e,uo,^'^), s and 
Uq G Be. 

Proof. As usual, we only give the formal derivation of this estimate which can be justified 
by approximating the variational solutions ui{t) and U2{t) by appropriate solutions of the 
regular equation (12. 9p . Set v{t) := Ui(t) — U2{t). Then, this function (formally) solves 

dtv = -A,{A,v - l{t)v + \v), dn{A,v - l{t)v + \v)\^ = 0, 
dfV — Arv + dnV + m{t)v = 0, on F, 



(5.16) 
where 



Kt)-= [ f'{sui{t) + {1 - s)u2{t))ds, m{t):= [ g' {sui{t) + {1 - s)u2{t)) ds. 
Jo Jo 

Multiplying this equation by {—Ax)^^v{t), integrating over Q and using the fact that 

{v(t)) = 0, we obtain 

(5.17) + W'^Mmhn) + mvit),vit))n < \Ht)\\l,^^^ + K\\vit)\\l,^^^, 
where K = ||fi''||c([-i,i])- We estimate the most comphcated term {l{t)v,v) as follows: 

(5.18) I l{t,x)\v{x)\'^dx> l_ l{t,x)\v{x)\^dx> 
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Thus, inequality fl5.17p reads 

(5.19) ~\\vmi^ + \\vMmiHn) + i^-mvmiHn)<^^^^^ 

Furthermore, using the trace-interpolation estimate 

\\v\\l2^r) < C\\v\\min)\\v\\LHn) < C(A - Xy'^HlHmin) + (A " >^)\Mhin)) 
and fixing 6 in such a way that KC < 1/2 (A (5) — A)^/^, we finally end up with 

(5.20) ||k(t)|||. + \\VMt)\\Un)+mvmh^n) + Mt)\\h^r)) < 2A||^T;(t)||i.(^), 

where /3 > 0. Using the Poincare inequality ||f ||j:/-i(f7) < C||f ||i2(Q), together with the 
Gronwall inequality, we deduce ( IS.lSp and finish the proof of the lemma. □ 

Thus, owing to Lemma [375| the semigroup S(t) is a contraction, up to the term \\9{ui — 
^2) ||l2([o,t]xQ)- The next lemma gives some kind of compactness for this term. 

Lemma 5.4. Let the above assumptions hold. Then, the following estimate holds: 

(5.21) \\dt{9iui - U2))\\L2([o,T],H-^n)) + 

+ \\9{ui - M2)||L2([0,T],ffi(n)) < Ce^^||Mi(0) - M2(0)||cI,™, 

where the constants C and K are independent of Ui{0) G B{s,Uo,^'^) , i = 1,2, and 
uo e Be. 

Proof. The second term in the left-hand side of fl5.2ip can be easily estimated by (15.71) 
(and the fact that VxO is uniformly bounded). So, we only need to estimate the time 
derivative. To this end, we recall that dtv {v = ui — U2) satisfies 

dtv = -Ax{A^v - l{t)v) 

in the sense of distributions. Therefore, for any test function (f G C^{Q), there holds 

{dtiOvit)), vp)^ = - (A,.^;(t) - lit)vit), A,.(M)f. = 

= (V.t;(t), V.A,(M)n + (KtHt), A,{9ip))^ . 

Since supp6' C ^5(^0), (15.111) yields 

I {l{t)v,Ax{9Lp))^ I < C\\v\\L2{n)\\^\\m{n) 

and, thus, 

I {dt{ev{t)),!f)^ I < Ci\\v{t)\\Hi{n)\\f\\H^n)- 

This estimate, together with (15.71) . give the desired estimate (15.211) on the time derivative 
and finish the proof of the lemma. □ 

It is now not difficult to finish the proof of the theorem. We introduce the functional 
spaces 

Hi := L'i[0,T],H\n))nH\[0,T],H-\n)), 

(5.22) 

^ ' m:= L\[0,T],L\Q)). 
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Then, obviously, Hi is compactly embedded into H. We also introduce, for any mq G 
the linear operator 

K„„ :B(£,Mo,$")^ei 

by 

K„„m(0) := eu{-), u{t) solves M 

(where the constants 6, T and the cut-off function 6 are such that Lemmas 15.31 and 15.41 
hold). Then, on the one hand, owing to Lemma [5.41 the map K^^ is uniformly Lipschitz 
continuous, 

(5.23) \\^uoiUi-U2)\\Mi < L\\Ul-U2h-", Ui,U2 e B{£,Uo,^'^), e<£o, 

where the Lipschitz constant L is independent of the choice of Uq G and e < Sq. On 
the other hand, it follows from Lemma 15.31 that 

(5.24) \\S{T)ui - S{T)u2U^ < (1 - 7)lki - + C\\K{ui - U2)\\m, 

where 7 > and C > are also independent of uq G Be, e < Eq and ui,U2 G B{e, uq, $"'). 

It is known (see, e.g., [I6j; see also |33j) that inequalities (15.231) and (15.241) . together 
with the compactness of the embedding Hi C H, guarantee the existence of an exponential 
attractor A4d{c) C B^ for the discrete semigroup S{nT) acting on the phase space Be 
(endowed with the topology of $"'). Furthermore, (15.71) . together with the control (15.51) 
of the time derivative, yield that the semigroup S{t) is uniformly Holder continuous with 
respect to time and space in [0,T] x Be. Thus, the desired exponential attractor A^(c) 
for the continuous semigroup S{t) on Be can be obtained by the standard formula 

M{c) := Ute[o,T]Md{c). 

Finally, although we have formally constructed the exponential attractor Ai{c) C Be C 
C°(f2) X C°(F) in the topology of only, the control of the C"-norm of Be, together 
with a proper interpolation inequality, give the finite-dimensionality and the exponential 
attraction in the initial topology of $c as well. This finishes the proof of Theorem l5.2[ □ 



6. Appendix 1. Some auxiliary results 

In this section, we establish several estimates which are used in the paper. We start 
with regularity results for the following singular elliptic boundary value problem: 

\A^u-u-f{u) = hi, in 
[ dn.u + u — Aru = h2^ on F, 

where hi G L'^{Q), h2 G -Z^^(F) and the nonlinearity / satisfies conditions (12.41) . As above, 
a solution u of this problem should be understood as a variational solution, analogously 
to Definition 13.11 Therefore, in order to justify the estimates given below, we factually 
need to deduce the corresponding uniform estimates for regularized problems of the form 
(16. ip . where / is replaced by its approximations /at (defined by (12. 7p ). and then pass to 
the limit 00. Since this passage to the limit is explained in details in Section [3], we 

give below the formal derivation of these estimates directly for the limit singular problem 
(16. ip . leaving the justifications to the reader. 
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Theorem 6.1. Let the above assumptions hold. Then, the following estimate holds for 
the solution u of problem f l6.1l) ; 

(6.2) II 

'"ll//i(n) + ll'"llHi{r) + Il^2||i2(r)), 

where the constant C is independent of hi and Furthermore, u G C°'{VL) nif^(r) with 
a < 1/4 and the following estimate holds: 

(6-3) Il«llc"(f7) + h\?H\T) < C{1 + \\hi\\l2(^n) + Il^2||i2(r)). 

Finally, F{u) E L\T), where F{z) := f{s) ds, S/^D^u eL'^{n),ue H'^{VL,), for every 
e > 0, where fi^ := {x G fi, f2) > e}, and the following estimate holds: 

(6.4) ||F(M)|Ui(r) + \\u\\]j2f^^^) + \\S/^Dru\\ + Il^2||i2(r)). 

Proof. Estimate (16.21) can be obtained by multiplying the first equation by u and inte- 
grating over VL (note that the existence and uniqueness of the solution u can be obtained 
exactly as in Section [3]). So, we only need to give a formal derivation of estimates (16.31) 
and n . 

The derivation of these estimates is based on a standard localization technique. Thus, 
we only give below a sketch of the proof, leaving the details to the reader. Let ^ be a 
smooth nonnegative cut-off function such that 0{x) = 1 if d{x, T) > e and 6{x) = if 
d{x, r) < e/2 which satisfies, in addition, the inequality 

(6.5) \v^e{x)\<ce^'^{x). 

Then, multiplying equation (16. ip by 

3 

^d^^{e{x)d^M), 

i=l 

integrating by parts and using estimate (16.21) (in order to estimate the lower-order terms) 
and the fact that /' > 0, we deduce that 

(6-6) \W\\]i\n,) < C{1 + \\hi\\l2in) + Il^2||i2(r)), 

where the constant C = depends on £ > 0, but is independent of u, hi and h2. 

Since C C°, a < 1/2, there only remains, in order to finish the proof of the theorem, 
to study the function m in a small e- neighborhood of the boundary F. 

Let Xo G r and y = y{x) be local coordinates in the neighborhood of Xq such that 
y{xo) = and Q is defined in these coordinates by the condition yi > 0. Then, in the 
variable y, problem (16. ip reads 

|^L=i ^y^(^ij(y)^y,^) + ELi Hy)dy,u + c{y)u - f{u) = hi, yi > 0, 
\EL=2 9yMj{y)9yjU) + ^J^2 e^Oy^u + g{y)u + h2 = dy,u, yi = 0, 

where aij, hi, c, dij, Ci and g are smooth functions which satisfy uniform ellipticity as- 
sumptions. 

Differentiating the first equation of (16. 7p with respect to yk, k = 2,3, multiplying the 
resulting equation by (pVk, where Vk := dy^u and is a smooth nonnegative cut-off function 
which is equal to one in the ball \y\ < e and zero outside the ball \y\ > 2e and satisfies 
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fl6.5p . using again the fact that /' > and the eUipticity assumption on the atj, we find, 
after standard transformations, 

(6.8) 7(0|Va;ffc|, |V^ffc|)n + i(t)Vk,aii{y)dy^Vk)r + < Ci\\u\\m{n) + Il^illi2(n))> 

where the positive constants C and 7 are independent of u. Differentiating then the second 
equation of ( 16 .7^ with respect to yk, inserting the expression for dy^^Vk thus obtained into 
fl6.8p and arguing analogously, we have 

(6.9) 'y{(j)\VxVk\, |Vxffc|)n + 7(0|Vy2,y3^fc|, \'^y2,ys'^k\)r + {4>Vk,Vk)n + i(pVk,Vk)r < 

< '^(IkllHi(n) + ll^illL2(n) + h\\m(r) + Il^2||i2(r))- 
Combining this estimate with (16.21) . we finally end up with 

(6.10) \\Ml^^^^,HHn,,y,)) + ll'^'"lli^«„^/HMg,.,3)) + \\My,=o\\Wl,.y,) ^ 

<C(l + ||/ii||i.(^) + ||/i2||i2(r)), 

where the constant C is independent of Xq and u. 

Returning to the variable x and using the fact that the boundary point Xq is arbitrary 
(and that F is smooth), we infer from (16.101) that 

(6.11) \\^xDTu\\l2(^a) + lkllH2(r) < C{1 + \\hi\\l2(^n) + Il^2||i2(r)). 
In addition, owing to the embedding 

L\R, H^{R'^)) n H\R, H^R"^)) C C''(R^), a < 1/4, 

estimate (16.101) . together with (16. 6p . also imply the estimate 

Ikllc-(Q) < C{1 + ||/ii||i2(f^) + WhWh^r))^ " < 1/4- 

Thus, in order to finish the proof of the theorem, we only need to estimate the L^-norm 
of F{u) on the boundary. To this end, we also use the localized equations (16.70 . but now 
multiply the first one by (pdy-^^u. Then, after obvious transformations, we have 

(6.12) dy,i^<l){y)au{y)\dyM' + 0(l/)^(w)) > 

> -C(0 + I V,0| + |D>|)(|V,«p + \dly^u\' + \dly.^u\' + \Dl,y,u\' + F{u) + l/^iH, 

where the constant C is independent of xq G F and u. Integrating this estimate with 
respect to ?/ G x R? and using (16. 2p and (I6.10p . together with the fact that 7^ only 
in a small neighborhood of the boundary, we see that 

(6.13) / (0(O,yi,?/2)i^(M(O,yi,l/2)) - 0(0, yi, 1/2)011(0, ?/i,?/2)|5yiM(0,?/i,?/2)nC?l/lC^Z/2 < 

< C{1 + ||/ii||i2(Q) + ||^2||i2(r)). 

Thus, keeping in mind the fact that (f){y) is nonnegative and is equal to one close to the 
given point Xq € F, we conclude, returning to the variable x, that 

\\F{u)\\L^(r) < C\\dnu\\l2^r^ + C(l + ||/ii||i2(^) + WhWh^r))- 
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There now only remains to note that, owing to estimate (16. lip and the second equation 
of (16. ip . we can control the L^-norm of (9„m on the boundary, 

ll^nM||i2(r) < C(l + ||/;.i||i2(f^) + ||/''2||L2(r))- 

The control of the L^-norm of F{u) on the boundary thus follows and Theorem 16.11 is 
proved. □ 

Our next task is to give an example of equations (16. ip for which the solution u does 
not satisfy the equations in the usual sense, but only in the variational sense described in 
Section[31 We recall that such an example cannot be found if the potential F{u) is singular 
at ±1. However, the situation is essentially different if the potential F{u) has finite limits 
as M — i> ±1. Indeed, in that case, the control of the L^-norm of F{u) is of no use and 
the singular part of the boundary (where \u{x)\ = 1) may now have positive measure and 
may even coincide with the whole boundary. As we can see from the following example, 
the equality [dnu]int = [dnu]ext can be violated at such singular points. 

Example 6.2. We consider the following example of a one dimensional boundary value 
problem of the form (16. ip : 

(6.14) y"-f(y) = 0, y'{±l) = K>0, a:G[-l,l], 

where the function / satisfies assumptions (12.40 and, in addition, -F(l) = Fi < oo and 
f{—y) = —f{y), which is of course a particular case of our general theory. Then, an 
analysis of the above ODE shows that, for relatively small values of K, this problem has 
a regular usual solution yK{x) which is odd, 

(owing to the symmetry and the uniqueness), and is separated from the singularities of 
/. However, there exists a critical value such that, for K > K+, yx coincides with 
the singular solution ?/+ of the problem 

yl-f{y+) = 0, y+{l) = l, y+{-l) = -1. 

Thus, the usual solution of (16.140 does not exist for K > K^. However, for these values 
of K, it can be uniquely defined as a variational solution. For the reader's convenience, 
we also give below a simple alternative proof of the above nonexistence fact which can 
be partially extended to the mult i- dimensional case. Since < 1, the usual 

interior regularity techniques (see the interior regularity estimate in Theorem 16.11) show 
that 

(6.15) \y'K{x)\<C, \yK{x)\<l-6, x G (-1/2, 1/2), 

where the positive constants C and 6 are independent of K. Multiplying now equation 
(I6.14P by y', integrating over [0, 1] and using (I6.15p . we obtain 

(6.16) \l\y'M\' - FiyKim < C, 

where the constant C is again independent of K. Thus, yx cannot satisfy the boundary 
condition y'xi^) = K if K is large enough and F{1) is finite. 
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Remark 6.3. Let yxix) be a variational solution of problem fl6.14p as constructed in the 
previous section. Then, since this solution is odd, it automatically satisfies the equation 

y"-f{y) = {y"-f{y)\-^,^] 

and, therefore, it is a (variational) equilibrium for the corresponding ID Cahn-Hilliard 
problem of the form fl2.ip . Thus, even in the ID case, problem (12.11) can have variational 
solutions which do not satisfy the boundary conditions in the usual sense. 
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